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We describe the creation of a Greenberger-Horne-Zeilinger 
(GHZ) state of the form (jOOO) + |111))/a/2 (three maximally 
entangled quantum bits) using Nuclear Magnetic Resonance 
(NMR). We have successfully carried out the experiment us- 
ing the proton and carbon spins of trichloroethylene, and con- 
firmed the result using state tomography. We have thus ex- 
tended the space of entangled quantum states explored sys- 
tematically to three quantum bits, an essential step for quan- 
tum computation. 

PACS numbers: 03.65.Bz, 89.70.+c,89.80.th,02.70.-c 

We live in a world which, to the best of our knowl- 
edge, is surprisingly well described by the laws of quan- 
tum mechanics. Fundamental processes in nature are all 
compatible with quantum mechanics, even if some of its 
predictions are counterintuitive. 

A good example of this peculiar behavior is given by 
the now famous pairs of two state systems described by 
Einstein, Podolosky and Rosen (EPR) Q . The unortho- 
dox behavior of this composite system has been crystal- 
lized by the Bell inequalities [||, which give a statisti- 
cal test for the existence of elements of reality ||. In 
our world, photons in a polarization state of the form 
(| 00) — \ll))/y/2 violate these inequalities as shown by 
Aspect et al. JiJ and therefore contradict the existence of 
elements of reality. 

In a beautiful paper, Greenberger, Home and Zeilinger 
(GHZ) H demonstrated that it is possible, in a single run 
of experiments using a state of the form 1 000) + 1 1 1 1) (the 
GHZ state), to refute the existence of elements of reality. 
The GHZ experiment has been succinctly summarized by 
Mermin [|| . Even though this experiment is very appeal- 
ing, nobody has yet been able to perform it. The reason 
is that it is rather difficult to precisely manipulate en- 
tangled quantum states of many particles. In fact, up 
to now, entangled pure states of only two particles have 
been systematically explored f|,@,||,D ■ 

In this letter we show how three particle entangled 
states can be realized using nuclear magnetic resonance 
techniques. We have carried out the experiment and ver- 
ified that indeed we had a GHZ state by using tomogra- 
phy pq |. We first describe how it is possible to obtain a 
pure state result from the initial mixed density matrix of 
a NMR system. Then we explain the sequence of opera- 
tions used to obtain a GHZ state and finally we give the 



experimental results. 

We will not investigate the non-local behavior of GHZ 
states, as NMR is not appropriate for this undertaking. 
What we have created, however, is a typical state needed 
for a three-bit quantum computer. Indeed our approach 
is the one used for quantum computation and we refer 
the reader to jll] for an introduction. 

The usual approach for quantum computation requires 
an initial pure state |Tl| ] . The computation itself is a 
transformation of this initial state using unitary oper- 
ations. Useful information is then extracted by a mea- 
surement of the final state. However, it has recently been 
shown that the computation can be performed using an 
initial mixed state as long as the decoherence time is suf- 
ficiently long @@0Jl|. 

In liquid state NMR, the computation takes place on a 
large ensemble of identical quantum systems. Each mem- 
ber of the ensemble of quantum systems consists of the 
interacting nuclear spins of a molecule in a high magnetic 
field. The initial state of the nuclear spins is achieved 
by allowing the system to relax to thermal equilibrium. 
Information processing with such an ensemble can be di- 
vided into three steps consisting of preparation, compu- 
tation and readout. Each of these steps is equivalent to 
an application of certain quantum operations identically 
to each member of the ensemble. 

The nuclear spins are manipulated by applying radio 
frequency (RF) pulses tuned to the Larmor frequencies 
of the spins 1611. The spins can be selectively excited by 
exploiting differences in Larmor frequencies. Entangle- 
ment and two spin operations are achieved by a delay to 
allow for interaction between spins. In our case, these 
interactions are scalar couplings which can be selectively 
turned off by the use of refocusing pulses tuned to one 
of the nuclei. Any unitary quantum operation can be 
decomposed into such operations [ p7| , ^8|Jl9| ] . 

The measurement step in NMR consists of observing 
the signal induced in RF coils by the precession of the 
nuclear spins |T^ ]. In effect, we measure the time evo- 
lution of the expectation of the o~ x and o~ y operators for 
each nuclear spin. By Fourier transforming the signal 
and analyzing the spectrum, other operators such as ten- 
sor products of either a x or a y with / or a z for pairs 
of interacting spins can be measured. The traceless part 
of the density matrix (called the deviation matrix) can 
be determined by a tomography procedure fici[| . This 
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involves several measurements of spectra after applying 
different reading pulses to the final state, thus permitting 
observation of all elements of the deviation matrix. 

The main problem in using NMR for the observation of 
multiparticle entanglements and quantum computation 
is that the sample is initially in a highly mixed state. 
There are several methods for overcoming this limitation 
without cooling the sample. One idea is to transform the 
initial mixed state so that we have a pseudo-pure state 

p pp =p\0...0)(0...0\ + ^-I, (1) 

where N is 2™ for n qubits and p+(l— p)/N is the prob- 
ability of the ground state (|0. . .)). In NMR this state 
is indistinguishable from a pure state because all observ- 
ables are traceless. The methods discussed in the litera- 
ture (T^Jl3|,[l4|,[l5| are all based on the idea of viewing the 
non-ground state components of the initial thermal den- 
sity matrix as noise and applying an averaging method 
to eliminate their contribution to measurements. These 
methods are in principle sufficient for exploiting the ad- 
vantages of quantum computing in an ensemble setting. 
In our experiments we used a new method for extracting 
pseudo-pure states from experimental data. The method 
is based on two features of our experiment. The first is 
that the sequence of pulses and delays applied to the sam- 
ple for preparing the GHZ state are much shorter than 
the decoherence time for the nuclear spins (about 5ms 
versus at least 210ms). As a result, the evolution is very 
close to unitary and therefore preserves the eigenvalue 
structure of the initial state. This structure is well known 
and determined by the thermal distribution. The unitar- 
ity property can be tested by comparing the eigenvalues 
of the output state to those of the thermal distribution. 
The second feature is that complete state tomography 
is used to analyze the output of the experiment. This 
allows us to decompose the output state into its eigen- 
states. The state |000) is associated with the smallest 
eigenvalue of the thermal state and must be transformed 
by unitary evolution into the corresponding eigenstate of 
the output. By comparing this eigenstate to the desired 
GHZ state, we learn whether our pulses indeed generated 
this state from that component of the ensemble which 
was initially in |000). 

In order to create a GHZ state we need three nuclei. 
A convenient system is trichloroethylene, the molecule 
shown in Figure 1. The spins of the hydrogen and car- 
bon nuclei were used for the three quantum bits. They 
interact only weakly with the chlorine nuclei, which can 
therefore be ignored. In a strong static magnetic field 
(11.5 Tesla along the z-axis) the evolution of the hydro- 
gen and the two carbon nuclei are well described by the 
Hamiltonian 

H = -who? - Wdcf 1 - uc 2 o° 2 + Jrc^g ^ + 
Jc lC2 (<xf 1 af 2 + af 1 af 2 + <rf 1 erf * ) + Jhc 2 of af 2 (2) 



with oj h w 500.1334915 MHz, u Cl ~ 125.7725805 MHz, 
loc 2 ~ 125.7732305 MHz, which gives a chemical shift of 
650 Hz between the carbons. The J couplings have values 
of J H d ~ 203 Hz, J Cl c 2 ~ 102 Hz and J HC . 2 w 10 Hz. 

The data was acquired with a Bruker DRX-500 
spectrometer using doubly labeled trichloroethylene 
( 13 C*i, 13 C* 2 ,99%). The relaxation time (7\) for the hy- 
drogen is 7s and 30s for the carbons. The phase deco- 
herence time (Ta) for hydrogen is 3s and 0.4s and 0.2s 
for G\ and Ci respectively. 

A simple circuit to create a GHZ state consists of 
a rotation by ir/2 around the y-axis followed by two 
CONTROL-NOT gates on the other qubits. Adapting 
this to our system results in the pulse sequence shown in 
figure 2 b). The deviation density matrix, initially given 
by the state 

p\ = uj H a h z (g) 1 ® 1 + u Cl 1 <8> erf 1 <8> 1 + ujc 2 1 ® 1 ® erf 2 (3) 

is transformed to 

p f A = -uJ H o h z <8> erf 1 <g> 1 - lu Ci o-^ <g> erf 1 <g> 1 

-wcfelSof 1 ®*? 3 (4) 

The tomography procedure was implemented with 
twelve sets of reading pulses. For each set, two exper- 
iments were performed to read off the hydrogen and the 
carbon spectra. 

To compute the peak intensities for each spectrum, 
we assumed that each peak is approximately Lorentzian. 
The peak positions and decay parameters associated with 
the Lorentzian shape were obtained by optimizing (in the 
least squares sense) the match to the calibration spectra. 
The carbon and hydrogen spectra were both matched to 
within 2%. This should be compared to the estimated 
noise, which was determined to be less than .5% for the 
carbon and .01% for the hydrogen nuclei. The mismatch 
is due to additional peaks, primarily from unlabeled com- 
pound, and also to shimming problems causing deviation 
from the ideal Lorentzian lineshape. The decay param- 
eters (essentially T 2 *) obtained after optimization were 
0.51 ± 0.03sec for the hydrogen nucleus, 0.41 ± O.Olsec 
for carbon 1 and 0.23 ± O.Olsec for carbon 2 (the error 
bars are estimated from the variation between different 
peaks of the same nucleus). The latter agree well with 
the experimentally determined T 2 time for the carbons. 
Since T 2 for the hydrogen is near 3sec, it can be seen 
that substantial peak broadening due to magnetic field 
inhomogeneity reduces the hydrogen T 2 . 

After the peak positions and decay parameters were 
determined, each spectrum of the experiment was decon- 
volved as a linear combination of the ideal peaks. The co- 
efficients of the combination yield the intensity and phase 
of each peak in the spectrum. These numbers were then 
normalized by the calibration intensities and phase cor- 
rected using the calibration phases. Then they were used 
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to determine the deviation density matrix of the output 
state of the experiment. 

An idea of how unitary the evolution has been 
is obtained by comparing the eigenvalues of the ini- 
tial deviation matrix ({24, 16, 16, 8, -8, -16,-16, -24}) 
to the final one ({24.1,16.3, 15.5,7.7, -8.1,-15.0, 
— 16.5,-23.8}). Clearly the transformation was unitary 
to a very good approximation . 

Using the procedure explained above we get, after di- 
agonalizing and taking the appropriate eigenvector, the 
experimentally determined GHZ state 



10" 



493 

-53 
-47 + 18z 
-28 + 5i 
25 + 22i 
-14- 31i 
-31 + 13i 



-53 
6 

5 - 2i 
3 

-3 - 2i 
2 + 3i 
3-i 



-47 - 18i 
5 + 2i 
5 
3 

-2 - 3i 
3i 
3 



-28 - 5z 
3 
3 
2 

-1 - 2i 



25 - 22i 


-14 + 3L 


-31 - 13i 


468 - 147i 


-3 + 2i 


2 - 3i 


3 + 2 


-51 + 16i 


-2 + 3i 


-3i 


3 


-39 + 3L 


-l+2i 


-2i 


2 


-25 + 13i 


2 


-2 


-2i 


30 + 14i 


-2 


2 


2i 


-23 - 25i 


2i 


-2i 


2 


-25 + 2L 


30 - Ui 


-23 + 25i 


-25 - 21i 


488 



A pictorial representation is given in figure 3. The fi- 
delity of the state compared to the ideal GHZ state is 



T= {^>GHZ\P 



GHZ 



GHZ) 



0.95. 



(6) 



In conclusion, we have shown how to construct an ef- 
fective GHZ state with a fidelity of 95% in NMR starting 
with the thermal mixed state. The experiments demon- 
strate the ability in NMR to fully explore the state space 
of multi-particle systems, which is all that is required 
for quantum computation. By using a four spin system, 
the paradoxical output of the originally proposed GHZ 
experiments can be observed as proposed in pOjj . How- 
ever, due to the microscopic separation of the particles 
involved and the method for observation used, this would 
not be a true test of the existence of elements of reality. 

To construct our GHZ state we developed a new 
method for extracting pseudo-pure states from NMR 
spectra. The method can be used efficiently for process 
tomography JlO| , since each experiment determines the 
transformation of the state space induced by the applied 
operation on each of the eigenstates of the input state. 

Our experiments demonstrate that room temperature 
liquid NMR is well suited for quantum computations in- 
volving small numbers of qubits. Although manipulating 
three qubits is a small step for large scale quantum com- 
putation, it is the first time that a quantum network 
has been used to systematically entangle more than two 
qubits. 
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FIG. 1. Pictorial representation of trichloroethylene. The 
three qubits are given by the nuclei of hydrogen (H) and of the 
two carbon 13's (Ci and C2). The last two are distinguishable 
because of the assymmetry of the chlorine environment. The 
chlorine nuclei have spin 3/2 and their interaction with the 
qubits can be neglected. 
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FIG. 2. a) Circuit to create a GHZ state. R a {9) corre- 
sponds to a rotation of the qubit by an angle 8 around the 
a-axis and the two other gates are CONTROL-NOTs (see Di- 
Vincenzo[14]). b) Implementation by an NMR pulse sequence 
producing a GHZ state (in the rotating frame of H and Ci). 
The delay di serves to let the coupling between the nuclei 
create the engtanglement, and (fe is used to generate a phase 
shift on the C2 nucleus. Coupling between the carbon nuclei 
during this delay is negligible. All pulses are non-selective. 
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FIG. 3. Pictorial representations of the theoretical and ex- 
perimental density matrices p for the GHZ state jOOO) + 1111). 
a), b) Real and imaginary parts of the theoretically deter- 
mined density matrix, c), d) Real and imaginary parts of the 
experimentally measured density matrix. Each bar graph rep- 
resents the value of the transitions between the eight states 
1 000} . . . I 111) by the height of the bar in the corresponding 
position of the 8x8 array. 
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